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ABSTRACT 

The characteristic of the solar acoustic spectrum is such that mode lifetimes get shorter and 
spatial leaks get closer in frequency as the degree of a mode increases for a given order. A direct 
consequence of this property is that individual p- modes are only resolved at low and intermediate 
degrees, and that at high degrees, individual modes blend into ridges. Once modes have blended 
into ridges, the power distribution of the ridge defines the ridge central frequency and it will 
mask the true underlying mode frequency. An accurate model of the amplitude of the peaks that 
contribute to the ridge power distribution is needed to recover the underlying mode frequency 
from fitting the ridge. 

We present the results of fitting high degree power ridges (up to £ = 900) computed from 
several two to three-month-long time-series of full-disk observations taken with the Michelson 
Doppler Imager (MDI) on-board the Solar and Hcliospheric Observatory between 1996 and 1999. 

We also present a detailed discussion of the modeling of the ridge power distribution, and 
the contribution of the various observational and instrumental effects on the spatial leakage, in 
the context of the MDI instrument. We have constructed a physically motivated model (rather 
than some ad hoc correction scheme) resulting in a methodology that can produce an unbiased 
determination of high-degree modes, once the instrumental characteristics are well understood. 

Finally, we present changes in high degree mode parameters with epoch and thus solar activity 
level and discuss their significance. 



Subject headings: Sun: oscillations 

1. Introduction 

Since from a single vantage point we can only 
observe a bit less than half of the solar surface, 
hclioseismic power spectra computed for a specific 
target mode with degree £ and azimuthal order m 
also contains power from modes with different - 
and usually nearby - £ and m values. The presence 
of these unwanted modes, or spatial leaks, compli- 
cates the fitting of the resulting observed spectra 



and degrade the mode parameter estimates, espe- 
cially when the leaks have frequencies similar to 
that of the target mode. As mode lifetimes get 
shorter and spatial leaks get closer in frequency 
(i.e., dv/d£ becomes small), individual p- modes 
can no longer be resolved. This mode blending oc- 
curs around £ — 150 for p- modes with frequency 
near 3.3 mHz and around £ = 250 for the f-modes. 
Once individual modes blend into ridges, as illus- 
trated in Figure 1, the power distribution of the 
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ridge masks the true underlying mode frequency 
and the ridge central frequency is not a good esti- 
mate of the target mode frequency. Moreover the 
amplitudes of the spatial leaks have been shown to 
be asymmetric (see for example Korzennik 1999). 
A direct consequence of this leakage asymmetry is 
to offset the power distribution of the ridges. Such 
offset results in a significant difference between the 
central frequency of the ridge and the frequency of 
the targeted individual mode. To recover the un- 
derlying mode frequency from fitting the ridge, an 
accurate model of the amplitude of the peaks that 
contribute to the ridge power distribution (i.e., the 
leakage matrix) is needed. 

The lack of unbiased determinations of mode 
frequencies and frequency splittings at high de- 
grees has so far limited the use of such high-degree 
data in helioseismic inversions for constraining the 
near-surface structure and dynamics of the sun. 
Since experiments like the Michelson Doppler Im- 
ager (MDI) on board the Solar and Hcliospheric 
Observatory (SOHO), with a two-arcsec-per-pixel 
spatial resolution in full-disk mode, allows us to 
detect oscillation modes up to £ « 1500 (Scherrer 
et al. 1995), only a small fraction of the observed 
modes are currently used. 

High-degree modes are trapped near the solar 
surface: for example, the lower turning point of a 
mode of degree I — 500 and frequency around 3 
mHz is at 0.99 of the solar radius. This makes 
them exceptional diagnostic tools to probe the 
near-surface region of the Sun, a region of great 
interest. Indeed, it is there that the effects of 
the equation of state are felt most strongly, and 
that dynamical effects of convection and processes 
that excite and damp the solar oscillations are pre- 
dominantly concentrated. Rabello-Soares et al. 
(2000) have shown that the inclusion of high- 
degree modes (i.e., £ up to 1000) has the poten- 
tial to improve dramatically the inference of the 
sound speed in the outermost 2 to 3% of the solar 
radius. Furthermore, inversion of artificial mode 
frequency differences resulting from models com- 
puted with two different equations of state (i.e., 
MHD and OPAL) recovered the intrinsic difference 
in Ti, the adiabatic exponent, throughout the sec- 
ond helium ionization zone and well into the first 
helium and hydrogen ionization zones, with error 
bars far smaller than the differences resulting from 
using two different equations of state. These tests 



were carried out using the relatively large obser- 
vational uncertainties resulting from ridge fitting, 
but with the implicit assumption that systematic 
errors were not present. These tests show that 
we can probe subtle effects in the thermodynamic 
properties of this region, but only when including 
such high-degree modes. 

The first estimates of high-degree mode fre- 
quencies used m-averaged Big Bear Solar Obser- 
vatory data (Libbrccht & Kaufman 1988). To re- 
cover the underlying mode frequency from fitting a 
given (n, £) ridge, they used a simple Gaussian pro- 
file as an approximation for the m-averaged leak- 
age matrix. Namely: 

C r 2 (£,£') = cx P (-(^_^) 2 ) (i) 

where A£ = £' — I while the e £ term represents the 
leakage asymmetry introduced by an image scale 
error of a fraction e of the image size. The ridge 
centroid frequency was estimated using a simple 
weighted average: 

. _ J2i'C?(£,£')A nr v nr 

where A n ^> is the individual mode power ampli- 
tude and v n ^i is the mode frequency. The fre- 
quency difference between the ridge and the mode 
frequency, Az/„ ;£ = v n j (ridge) - ^(mode), can 
thus be estimated using a parametric representa- 
tion of the leakage coefficients. 

For intermediate degree modes (50 < I < 150) 
this frequency difference can be directly measured 
by reducing the frequency resolution of the ob- 
served power spectra as to force individual modes 
to blend into ridges. Using observed Av, the pa- 
rameters s and e can be calibrated at these inter- 
mediate degrees and the correction extrapolated 
to high-degree modes. 

As Libbrecht & Kaufman (1988) state in their 
paper, this is only a first step and there is sub- 
stantial room for improvement particularly at high 
I. Following the same idea, but improving on 
the method, Korzennik (1990) and Rhodes et al. 
(1999) estimated high-degree mode frequencies us- 
ing Mount Wilson and MDI data respectively. In 
contrast, Bachmann et al. (1995), using data from 
the High-L Helioseismometer at Kitt Peak, calcu- 
lated the m-averaged leakage matrix, but neglect- 
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ing the horizontal components, in order to esti- 
mate high-degree mode frequencies. 

In this paper, we present the results of an exten- 
sive study of the various elements that contribute 
to the precise value of the effective leakage ma- 
trix, that in turn is key to the precise determina- 
tion of high degree modes. We have attempted to 
construct a physically motivated model — rather 
than an ad hoc correction scheme — in order to 
produce an unbiased determination of the high- 
degree modes. Since Korzennik (1999) has shown 
that the inclusion of the horizontal component of 
the leakage matrix calculation partially explains 
its observed asymmetry, we have included the hor- 
izontal component in all our leakage matrix calcu- 
lations. 

In Section 2, we describe the data we used and 
how we computed and fitted the power spectra 
used in this work. In Section 3 we present and dis- 
cuss ridge modeling for high degree modes while 
in Section 4 we address the issue of estimating 
the ratio between the radial and horizontal com- 
ponents. In Section 5 we discuss the instrumental 
effects specific to MDI that must be included to 
properly model the ridge power distribution. Fi- 
nally, in Section 6, we present the results of our 
analysis followed by our conclusions. 

2. Data Used and Fitting Methodologies 

The data used in this work consists of four sep- 
arate time series of full-disk Dopplergrams. These 
Dopplergrams were obtained by the MDI instru- 
ment while operating in its full disk 2" per pixel 
resolution mode and during the Dynamics Pro- 
gram. The instrument is operated in this mode 
some 3 months every year, when the available 
telemetry bandwidth is large enough to bring 
down full-disk images. Table 1 lists the starting 
time and duration of each of the four Dynamics 
data time series. 

The spherical harmonics decomposition up to 
t = 1000 was carried out by the MDI Science 
Team (Schou 1999). For the same periods, mode 
frequencies and rotational splitting coefficients of 
individual p-modes were computed for low and in- 
termediate degrees as part of the Structure Pro- 
gram (Schou 1999). 

For each Dynamics Program period we com- 
puted and fitted power spectra using two distinct 



approaches. The purpose of the first one was to 
get some observational estimate of the leakage ma- 
trix asymmetry, using an optimum estimate of the 
sectoral limit spectra computed with a high fre- 
quency resolution where individual p-modes can 
be clearly resolved. 

We used only sectoral spectra because the leak- 
age matrix for the sectoral modes is diagonal at a 
fixed radial order, n (namely non-zero only when 
5m = 5£), resulting in a simple leakage pattern, 
with modes separated by at least the quantity 
dvjdl. For zonal and tesseral modes the leakage 
pattern is more complex resulting in leaks very 
close to each other, some separated only by even 
multiples of dvjdm. For these modes, only for 
very low frequencies (where the lifetime is long 
enough), and if the time series is long enough, can 
the main peak be resolved from its closest spatial 
leaks. 

The second approach is aimed at computing op- 
timum power spectra for high-degree modes ridge 
fitting. Namely spectra with a lower frequency res- 
olution — since individual modes are not resolved 
- while in the process averaging out the realiza- 
tion noise that a high frequency resolution spec- 
trum would reveal. Moreover, at low and interme- 
diate degrees individual modes blend into ridges 
as a result of the low frequency resolution. This 
approach creates an overlap between mode fitting 
and ridge fitting where the frequency offset intro- 
duced by ridge fitting can be directly measured. 
By using this overlap, we can check the model we 
have developed to infer mode parameters from the 
fitted ridge parameters. 

2.1. Limit Spectra for Resolved Modes at 
Intermediate Degrees 

Estimates of limit spectra were computed for 
intermediate degree sectoral modes, when individ- 
ual modes could be resolved, using the following 
procedure. First, each time scries of spherical 
harmonic coefficients was detrended using a 21- 
minute-long running mean. Then, the 9th order 
sine multi-tapered power spectrum was computed, 
using the two- to three-month-long time series and 
thus corresponding to a frequency resolution be- 
tween 0.18 /iHz and 0.24 /aHz. Finally, for a range 
of degrees (15 consecutive £) at a fixed order, n, 
small fractions of power spectra (about 300 ^Hz), 
centered around each mode frequency, were aver- 
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Table 1 
Dynamics Program EPOCHS 



Data Set 
[year] 


Starting Date 
month/day 


Duration 
[days] 


1996 


05/23 


63 


1997 


04/14 


93 


1998 


01/09 


92 


1999 


03/13 


77 



aged. The sectoral mode frequencies were com- 
puted using a table of frequencies and frequency 
splittings resulting from mode fitting of MDI data. 

A seven-component profile was fitted in the 
least-squares sense to these limit spectra using the 
following parameterization: 



£+3 



1 + a n ,i(x - a nt t/2) 



£'=£-3 
+ b„,£ + S n j {v - U nt t) 



1 + x 2 



where 



V-Vr, 



ln,£ 



(3) 



(4) 



and where A n ^i and i/ n j> are the individual mode 
amplitude and frequency respectively; ^ n ,£ and 
a n ,e are the mode width and asymmetry param- 
eter respectively - assumed constant for all seven 
components; and b n ,£ and s n ,e represent the back- 
ground power. 

Note that the shape of the profiles used in this 
formulation is equivalent to the one defined by 
Equation 4 of Nigam & Kosovichev (1998). After 
some rudimentary algebra, one can identify their 
asymmetry coefficient, B, to the one we use, a, 
namely § = B {I - ^-). 

2.2. Power Spectra at High Degrees and 
Ridge Fitting 

For each period listed in Table 1, the time se- 
ries of spherical harmonic coefficients were sub- 
divided in 4096-minute-long intervals. For each 
spherical harmonic degree, £, and each azimuthal 
order, m, each interval was first detrended using 
a 17- minute-long running mean, and then Fourier 



transformed with an oversampling of 2 to estimate 
the power spectrum. The power spectra corre- 
sponding to each 4096-minute-long interval in a 
given period were then averaged together to pro- 
duce high signal-to-noise ratio spectra with a fre- 
quency resolution of 4.1 /iHz. As a result of this 
low frequency resolution, the individual modes 
down to I w 100 blended into ridges, while the 
ridge width itself is still very much oversampled. 

The following asymmetric profile plus a back- 
ground term was fitted in the least squares sense 
to each low frequency resolution spectrum, for 
100 < £ < 900: 



Pf„ m {v) 



where 



and 




1 + OLn,l,m{x n ,l,m ~ &n,t,m/2) 



l + X 



52 
n.i.ra 



(5) 



Xn : £,m 



v — Vn,e,m 
Jn,£,m 



log B t . m (v) = ^2 bj(£,m) v 3 

3=0 



(6) 



(7) 



such parameterization of the background power is 
better suited for the fit performed here, namely 
over a broad frequency range. 

In this case, A n ^ m , 7„^, TO , J>„,£, m and a n> t >m 
are the ridge amplitude, width, frequency and 
asymmetry parameters for a given (n, £, m) ridge 
fitting. For practical reasons, this fitting was car- 
ried out only every 5th £ for 100 < £ < 250 and 
every 10th £ for 250 < £ < 900, and only for 
some 50 equally spaced m values at each £. From 
the resulting ridge frequencies, the ridge frequency 
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splittings were parameterized in terms of Clebsh- 
Gordon coefficients (a,, for i = 1,6, see Ritzwoller 
& Lavely 1991). 

A single asymmetric function was used to to fit 
each ridge power density profile (Eq. 5). Since the 
ridge power density profile results from the over- 
lap of the target mode and the spatial leaks, it 
is mathematically not a similar single asymmet- 
ric function. In the specific case of high degree 
modes and at the frequency resolution we have 
used we believe that our approach is both jus- 
tified and practical, based on the following con- 
siderations. First, the profile resulting from the 
overlap of nearby profiles is reasonably well mod- 
eled by a single profile when the ratio of the width 
of the individual profiles (r e fj) to the separation 
(dv/dl) become large. Indeed, for a ratio of 2 - a 
ratio that roughly corresponds to the mean ratio 
around £ = 300 - the root-mean-squares (RMS) of 
the difference between a superposition of profiles 
and a single profile is only 4% (when normaliz- 
ing the resulting fitted profile maximum to unity). 
For a ratio of 4 (£ « 500) the RMS of this differ- 
ence drop to 2%, while for a ratio of 10 (£ « 700) 
this RMS drop to 0.6%. These numbers are to 
be compared to the RMS of the residuals to the 
fit that are around 5 to 7% (when using the same 
normalization) at all degrees and are dominated 
by the realization noise. The second reason to use 
a single profile is practical: fitting to the resulting 
power density profile the sum of individual over- 
lapping profiles is underconstrained, unless some 
assumptions on the separations and/or the rela- 
tive amplitudes are made. Since neither of these 
quantities are known a priory with enough preci- 
sion such approach is impractical. 

3. Ridge Modeling 

We describe here our physically motived ridge 
model, from which corrections to the ridge pa- 
rameters are estimated. Once the offsets between 
ridge and mode parameters are determined, the 
observed ridge parameters can be corrected and 
the mode parameters for high-€ obtained, making 
it possible to study the near-surface region of the 
Sun. 



3.1. Leakage Matrix: Radial and Horizon- 
tal Components 

Since the spatial dependence of solar oscilla- 
tions can be described in terms of spherical har- 
monics Y™(6, 4>) of co-latitude 9 and longitude 0, 
the velocity signal resulting from these oscillations 
can be written as: 

VnAm = (V r YP, V h dgY e m , V h C^Y/"). (8) 

and the leakage matrix components, resulting from 
the spherical harmonic decomposition of the ob- 
served line-of-sight velocity, are given by: 



r 9 



<j> WY e m *Yp' sin 9 cos cf> dfl (9) 
- j I WYl n *d Y^' cos 9 cos <f> d®0) 



WY^d^YJJ 1 smfidn (11) 



where * represents the complex conjugate opera- 
tor and W{9, 4>) is the spatial window function of 
the observations - i.e., a function that delimits 
the angular span of the observations and that in- 
cludes additional spatial attenuations like spatial 
apodization. 

The complete leakage matrix is the sum of the 
radial component C r and the horizontal compo- 
nents C e + C^. Thus a mode with a given (£' , to') 
will leak in the (£, to) power spectrum with an am- 
plitude attenuated by a factor Cj m . t , m , where: 



and where /3„ ; £ is the horizontal-to-vertical dis- 
placement ratio. Using a simple outer boundary 
condition, i.e. that the Lagrangian pressure per- 
turbation vanish (5p = 0), the small amplitude 
oscillations equations for the adiabatic and non- 
magnetic case lead to an estimate of the ratio (3, 
given by (Christensen-Dalsgaard 1997): 



G Mm L 



Rq 



n,e 



21 



(13) 



where G is the gravitational constant, M Q is the 
solar mass, Rq is the solar radius, u is the cyclic 
frequency (u) — 2ttv) and L 2 — £(£ + 1). It can be 
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reduced to the square of the ratio of the frequency 
of the fundamental (n = 0) for the given degree 
to the mode frequency, since the frequency of the 
fundamental is given by (Gough et al. 1980) 



2 _G M Q L 



(14) 



Since for a given £, the f-mode frequency is smaller 
than any p-mode frequency, Equation 13 implies 
that < f3 < 1. 

Estimates of (3 resulting from computing the 
adiabatic and non-magnetic cigcnfunctions using 
a standard solar model (Christensen-Dalsgaard 
et al. 1996) at the appropriate atmospheric height 
for the MDI instrument and using a different up- 
per boundary condition than the one leading to 
Equation 13 gives a very similar result. These 
calculations produce values that, above 2 mHz, 
decrease with frequency slightly faster than pre- 
dicted by Equation 13 (i.e., 15% smaller at 3 
mHz). However, the non-adiabatic nature of the 
oscillations in the solar atmosphere are not taken 
into account in these calculations nor is it taken 
into account in the derivation of Equation 13. 
Since the solar interior is non-adiabatic near the 
surface and since the effect of the magnetic field, 
if any, is likely to be mostly near the surface, the 
real value of (3 might be somewhat different. 

The estimate of the coefficients of the leakage 
matrix given by equations 9 to 11 is a useful guide 
but remains oversimplified. Indeed, these equa- 
tions do not take into account the pixel size of 
the detector or the instrumental modulation trans- 
fer function (MTF), effects that are substantial at 
high degrees. Therefore, rather than integrating 
these equations, we opted to compute the leak- 
age matrix coefficients by constructing simulated 
images corresponding to the linc-of-sight contri- 
bution of each component of a single spherical 
harmonic mode and decompose that image into 
spherical harmonic coefficients using the numer- 
ical decomposition used to process the observa- 
tions. Such an approach allows us to easily include 
effects like the finite pixel size of the detector as 
well as other instrumental effects like the instru- 
mental MTF, a plate scale error or some image 
distortion. 

Since this calculation is computationally ex- 
pensive and since the leakage coefficients vary 
smoothly with I and to, it was carried out only 



for a subset of modes, from which all the required 
values were interpolated. 

Figure 2 shows some typical examples of ob- 
served and synthetic spectra where the amplitude 
of the leaks corresponds to leakage matrices cal- 
culated with and without including the horizon- 
tal component (right and left panels respectively) . 
This figure clearly shows not only that the ob- 
served amplitudes of the spatial leaks are asym- 
metric, but the importance of the contribution of 
the horizontal component in the leakage matrix 
calculation, since the radial component of the leak- 
age matrix is nearly symmetric. It should also be 
noticed that, although the complete leakage ma- 
trix agrees in broad terms with the observed spec- 
tra, it does not agree in detail, as already pointed 
out by Korzennik (1999). This was the motivation 
for the work we present here, where we discuss in 
details how to best compute the correct effective 
leakage matrix, i.e. the one resulting from the ac- 
tual data analysis. 

3.2. Perturbation of the Leakage Matrix 
by the Solar Differential Rotation 

Woodard (1989) has computed the distortion 
of the eigenfunctions by a slow, axisymmctric dif- 
ferential rotation. He showed that the distorted 
eigenfunctions can be expressed as a superposition 
of the undistortcd ones. Under the assumption of 
axial symmetry, each of the superposed functions 
has the same azimuthal order to, while symmetry 
with respect to the solar equatorial plane implies 
that the values of the spherical harmonic degree 
included in the superposition be either all odd or 
even. The Coriolis acceleration is also assumed to 
be small which is valid for high-degree p-modes 
(£ > 100). This additional simplification reduces 
the perturbed eigenfunctions to be approximately 
a superposition of unperturbed eigenfunctions of 
the same radial order n. 

Also, the rotation rate for high degree modes 
(£ > 100 corresponds roughly to the radius range 
O.9i? < r < Rq) is primarily a function of lati- 
tude, and can be parameterized as: 

Q(6>) = 2tt(B {) + B 2 cos 2 9 + B t cos 4 6) (15) 

where Bq — 473 nHz, B 2 = —77 nHz and B4 = 
—57 nHz (Snodgrass & Ulrich 1990), and 6 repre- 
sents the co-latitude. 



6 



The perturbed leakage matrix can be expanded 
in terms of the unperturbed one as: 



(16) 



where: 



< l , = f cos [p0 + d9 (17) 



<f$(0) 



2tt 



1 



2 cV<9£ 
1 



B 4 



/ ) sin 9 



-^-B iy 4 e ,sm(29) 
lb 



(18) 



I" is even, 



and where Ggi e<< is zero unless 
2p = £'-£", yj = l-(m/L) 2 . 

This distortion is illustrated in Figures 3 and 
4, where we show the variation of Gci" with re- 
spect to m and I — I' for selected values of I, and 
a synthetic spectrum with and without distortion 
by differential rotation for I = 500. We should 
also add that Woodard (2000) has shown that in- 
cluding meridional circulation affects these calcu- 
lations at the 10% level. 

3.3. Theoretical Ridge Parameters 

Synthetic power spectra were computed by 
overlapping individual modes as they leak in the 
simulated spectrum according to: 



P n ,t,m{v) 



e+3 m+3 

^ X] Gi m . e , m ,(/3 n ,e) 

£'—£—3 m'=m— 3 
1 + 6i n ,l{x n ,l' ,m' ~ 6tn,t/2) 



1 + %n,l'm' 



v - v„ 



%n,£' ,m' 



(19) 



(20) 



where v, 7 and a are the mode frequency width 
and asymmetry parameter used to generate the 
model. We used for (z>, 7, a) values based on 
observed parameters and used Equation 13 for the 
values of j3. We generated these synthetic spectra 
for values of I between 100 and 900 - but only 
every 5th t up to 250 and every 10th t above 250, 
replicating the sampling we used in the fitting of 
the observed MDI spectra as described in Section 
2.2. 



We generated two sets of synthetic power spec- 
tra, one where we did not include the effect of the 
distortion of the eigenfunction due the differential 
rotation (i.e., Gi>.i" — Si'.i"), while the other did 
include this effect. We then used on these two 
sets of simulated ridges the same fitting procedure 
we used on the actual observations, as described 
in Section 2.2 and Equation 5, to compute ridge 
frequencies and frequency splitting coefficients. If 
our model is correct and complete, our simulations 
should produce frequency and splitting coefficient 
offsets smilar to the one resulting from MDI ob- 
servations. 

Figure 5 illustrates, in an l-v diagram format, 
the region of overlap between mode fitting and 
ridge fitting. This range of frequencies and de- 
grees is useful to valide our methodology but too 
restrictive to dare to carry out any extrapolation. 
In figures 6 to 8 we show frequency and split- 
ting coefficient offsets from our two sets of sim- 
ulations and the offsets from actual observations 
for intermediate-degree modes where there is an 
overlap between ridge and mode fitting. The the- 
oretical differences agree overall with the observa- 
tions (i.e., magnitude and general trend), but not 
in detail. The effect of the distortion by the differ- 
ential rotation is significant, and affects the most 
the odd splitting coefficients. It explains very well 
the discontinuity in the splittings coefficients seen 
as far back as in Korzennik (1990), and more re- 
cently in Rhodes et al. (1999). 

3.4. Influence of the Rotation Profile on 
the Splitting Corrections 

The latitudinal dependence of the solar rotation 
rate, £1(0), must in principle be known to compute 
the perturbation of the leakage matrix by the dif- 
ferential rotation (see Section 3.2). Equations 17 
and 18 result from an explicit parameterization of 
the near surface rotation with latitude given by 
Equation 15, and thus depend on the precise val- 
ues adopted for B 2 and B4. Therefore the cor- 
rection for the ridge rotational splittings depends 
on a model of rotation, while the corrected split- 
tings are meant to allow us to derive that rotation 
profile. 

In order to estimate the effect of the rotation 
profile on the correction of the splittings, we have 
synthesized theoretical ridges using three different 
rotation profiles. Each set of ridges was then fit- 
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ted and ridge frequency splittings were computed, 
from which three sets of splitting coefficient cor- 
rections were computed. 

The first model, hereafter model A, used a 
somewhat arbitrary rotation profile defined by 
B 2 = — 75 and B4 = —50 nHz, model B used a 
profile inferred by Schou et al. (1998) using the 
MDI medium-^ splitting coefficients, while model 
C used the profile given by Snodgrass & Ulrich 
(1990) based on the cross-correlation of Mt. Wil- 
son magnetograms. The three profiles are com- 
pared in Figure 9. Relative variations with respect 
to model A are show in Table 2 and illustrated 
in Figure 10. The first column of Table 2 indi- 
cates the average relative change in the rotation 
profile, while the other columns show the average 
and RMS of the respective parameter relative cor- 
rections. 

This table shows that the frequency correction 
is affected by the choice of the rotation profile, 
but at a relative level substantially smaller than 
the relative change in the rotation profile. We 
should also note that the relative change in the 
frequency correction is a function of degree, as it 
grows larger at larger degrees. Table 2 also shows 
that the correction in the cti coefficients is affected 
by the choice of by a relative amount compa- 
rable to the relative change in the rotation profile. 
Again we should point out that the change shows 
a variation with degree and frequency, ft is most 
pronounced for a\ for the f-mode, but is also sub- 
stantial for n = 1, while for 03 only the f- modes 
show a substantial effect. Of course some of this 
is the direct result of the variation of AO with lat- 
itude between the different models used. Finally 
we should point out that the solar rotational pro- 
file is known with a better precision than the num- 
bers listed in Table 2 suggest, resulting in smaller 
effects. 

4. Estimate of (3 from Intermediate Degree 
Observations 

Korzennik (1999) presented a methodology to 
estimate (3 from intermediate degree observations 
and speculated, on the basis of a single season of 
Dynamics observations, that Equation 13 might 
not be satisfied. This methodology proceeds as 
follow: First, the leakage asymmetry is quantified 



using a single factor, S x , defined by: 

S x (n,£ 



using the values of A n ^i fitted to the sectoral limit 
spectra as described in Section 2.1 and Equation 3. 
Next, a correspondence between S x and (3 is estab- 
lished using the complete leakage matrix, namely: 

S x (f3, I) = ^-^Y ' (22) 

Finally, using this correspondence (shown to be 
nearly independent of £ in Figure 4 of Korzennik 
1999) the values of S x are converted to observa- 
tional estimates of (3, hereafter refered as (3. 

Alternatively, one can estimate (3 using a for- 
ward approach, where f3 is adjusted as to produce 
theoretical ridge frequencies — by fitting synthetic 
spectra — that matches the observed ridge fre- 
quencies. Both methodologies agree, as shown in 
Figure 11, for a single year of Dynamics data. 

When this method was applied to successive 
years, using the same leakage matrix, we found 
that f3, while agreeing for 1996 and 1997, changed 
noticeably in 1998 and substantially in 1999. Since 
it is highly unlikely that (3, the horizontal-to- 
vertical displacement ratio, would actually change 
over the years, we came to the obvious conclusion 
that observed discrepencies — whether expressed 
in term of leakage asymmetry or frequency or split- 
tings differences — resulted from some missing 
component in our modeling. 

Thus instead of trying to adjust the leakage ma- 
trix with some ad hoc estimate of (3, one must 
improve the leakage matrix estimate using phys- 
ically motivated contributions, as to produce an 
unbiased determination of the high-degree mode 
frequencies. In fact, we show in the next sections 
that the differences of (3 between 1998 and 1999 re- 
sults from a time variation of the plate scale error 
in the decomposition of the MDI images. 

We should also point out that Schou & Bog- 
art (1998) while investigating flow and horizon- 
tal displacements from ring diagrams, using MDI 
Dynamics observations, found the ratio (3 to be 
in good agreement with the theory (see Figure 
7 of Schou & Bogart 1998). Schmidt, Stix, & 
W6hl (1999) have argued that the observed ve- 
locity vectors in the photosphere appeared to be 



(21) 
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Table 2: Relative variation of ridge to mode corrections with respect to the rotation profile. 





An 

[%] 


Ay 

[%] 


Aai 
[%] 


Aa 3 
[%] 


Aa 5 

[%} 


Model B 


- A 


—AA 


+0.22 ±0.05 


-2.52 + O.fO ^ 


+f.28±0.09 W 


+f0.55±0.08 


Model C 


- A 


+6.2 


-1.65 + 0.52 


+6.47 + 0.06 


+7.16 + 0.06 


+8.28 + 0.05 



a computed for n > 1 only 
b computed for n > only 



more vertical than predicted, using MDI velocity 
power spectra to study the centre-to-limb varia- 
tion. But Rhodes et al. (2001) concluded that the 
ratio /3 is very close to the theoretical value, by 
studying m-averaged intermediate degree GONG 
power spectra. 

5. Instrumental Effects 

The MDI instrument has been very stable over 
the years. Nonetheless, continuous exposure to so- 
lar radiation has increased the instrument front 
window absorption resulting in a long term in- 
crease in its temperature (Bush, Chu, & Kuhn 
2001). Moreover, there is indirect evidence for a 
temperature gradient from the center to the edge 
of the front window, associated with that temper- 
ature increase. This gradient is believed to cause 
a small curvature that converts the front window 
into a weak lens. 

Evidence of this is a drift in the instrument fo- 
cus (see top panel of Figure 12). The short term 
small variation in the focus is correlated with the 
annual temperature change of the front window 
due to the satellite orbit around the Sun (see Bush, 
Chu, & Kuhn 2001, for further details). 

Image focusing on the MDI instrument is con- 
troled by inserting (or not) blocks of glass of vary- 
ing thickness into the light path (see Scherrer et al. 
1995). There are 9 possible focus positions, with 
one focus step corresponding to approximately a 
third of a wave. Most of the time, the image was 
slightly defocused on purpose to match the detec- 
tor resolution. The actual configuration in which 
the instrument was operating is shown in the top 
panel of Figure 12, with the corresponding amount 
of defocus. 

A direct consequence of a change in focus is a 
variation of the image scale at the detector. The 



bottom panel of Figure 12 shows the variation of 
the ratio of the radius calculated using the ini- 
tial plate scale relative to the mean image radius. 
The mean image radius is the instantaneous aver- 
age of the image semi-major and semi-minor axis 
(see below), while the initial plate scale refers to 
a fixed estimate of the image radius computed in 
early 1996. The discontinuities arc due to changes 
in the instrument focus position and are indicated 
by vertical dashed lines. The observed annual vari- 
ations as well as the small but systematic increase 
with time are correlated with the aforementioned 
changes in the front window temperature. 

Besides the image scale, another relevant aspect 
of the optical characteristics of the instrument is 
the distortion of the solar image. The shape of the 
solar limb on the detector can be approximated 
by an ellipse with an eccentricity 100 times larger 
than the actual eccentricity of the sun. The dif- 
ference beween the semi-major and the semi-minor 
axis of the MDI full-disk image is 0.50 ± 0.04 pix- 
els, while the mean radius is 488.45 pixels. The 
orientation is such that the semi-major axis forms 
an angle of 34° ± 1° with the equator in the clock- 
wise direction. Although the changes in the MDI 
focus affect the solar image ellipticity, these varia- 
tions are rather small, resulting in a distorted but 
relatively stable image. 

Let us point out that the spherical harmonic 
decomposition of all the Dynamics data has been 
carried out assuming a circular solar limb, with 
no image distortion. The image radius used for 
this decomposition was computed using the value 
corresponding to the initial estimate of the MDI 
instrument plate scale and rescaled according to 
the distance to the sun given by the spacecraft or- 
bit ephemcris. The variations introduced by the 
change of focus position and the effect of the front 
window were not all known or well characterized 
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and thus were not taken into account at the time 
the spatial decomposition was computed. Since 
the computation of 500,000 coefficients for some 
468,000 images - each 1024x1024 pixels - is a 
rather intensive computing task, the spatial de- 
composition has not yet been recomputed to take 
into account the known image scale variations. 

5.1. Effects of Plate Scale Error 

Variations of the mean image radius relative to 
the initial plate scale, even at the 10 -3 level (as 
illustrated in the lower panel of Figure 12) will 
somewhat affect the amplitude of the spatial leaks 
(see Equation 1) and must be taken into account 
in the leakage matrix calculation. 

In Figure 13 we compare the computed limit 
spectrum (see Section 2.1) estimated for n — 2 and 
I = 103 for 1998 and 1999, in the upper panel. It 
clearly shows that the leakage asymmetry changes 
sign between both epochs (Ag-i < A? + i in 1998, 
while Az-x > A^ + i in 1999). The lower panel 
of Figure 13 shows simulated power spectra com- 
puted using a leakage matrix with and without a 
plate scale error of —0.1%. As expected 1 , the plate 
scale error causes a change in the leakage asymme- 
try similar to the change observed between 1998 
and 1999 observations. 

The plate scale errors for the different Dynam- 
ics epochs are given in Table 3. Note also that it 
is the size of the equatorial radius that will dictate 
the amplitude of the leaks for the sectoral spectra. 
Since the shape of the solar image on the detector 
is nearly an ellipse, the radius at the equator is 
slightly larger (« .02%) than the mean radius. 

The effect of the plate scale error can be clearly 
seen in the observed ridge frequency offsets pre- 
sented in Figure 14. The offsets between ridge and 
mode frequencies for 1996 and 1997 are very simi- 
lar, while for 1998 and especially for 1999 they be- 
come substantial. The values of Aai (i.e., the first 
ridge splitting coefficient offset) also show some 
variation from year to year, although smaller. 

We have modeled the effect of a plate scale 
error by computing the leakage matrix resulting 
from decomposing images using the wrong radius, 



and including the effect of the distortion of the 
eigenfunction by the differential rotation. These 
leaks were then used to generate synthetic power 
ridges that were fitted using the same ridge fit- 
ting procedure described earlier. Figure 15 shows 
the frequency offsets and rotational splittings first 
coefficient offsets for different plate scale errors. 
Such simulations reproduce rather well the varia- 
tions observed in the data, including the changes 
between epochs. Note that at this point we have 
not yet included any image distortion. 

The effect of a small plate scale error on the 
ridge frequency can be approximately modeled as 
a wavenumber scaling error in the spatial decom- 
position, namely £ € = £(l — e), causing an apparent 
frequency shift of the ridge. Using a linear approx- 
imation, valid for small e, this can be written as: 

Au = u(Q-u(£) = ^A£=-^e£ (23) 

In practice the observed frequency differences 
result from the combination of a plate scale error 
and some other terms. Thus, the variation of this 
frequency difference with epoch, if the only change 
is a plate scale error, can be measured, and corre- 
sponds to: 

d ^ f ^Sff-^Sff (24) 



Ai/(P 3 ) - Ai/(Pi) 
e(P 2 ) - eJA) 
e(P 3 )-e(P 1 ) 



(25) 



1 ln this case, the leakage asymmetry is the result of two fac- 
tors: the contribution of the velocity horizontal component 
that leaks more power at higher £ and a negative plate scale 
error that has the opposite effect. 



where e(P») is the plate scale error for the epoch 
Pi- 

Similarly, the variation of the rotational split- 
tings first coefficient differences, Aai, due to a 
plate scale error will be proportional to first or- 
der to e. Our simulations have shown a residual 
dependence on frequency and degree, namely: 

A ai (P 2 ) - Aai(Pi) <* A (e(P 2 ) - e(Pi)) (26) 

where 

A = Ao + A v {v - v ) + Ae(£- £ D ) (27) 

and A = —580 nHz, while A v = — 26 nHz/mHz 
and At = 0.13 nHz, for v Q = 3.333 mHz and £ = 
150, using our normalization of the Clebsh-Gordon 
coefficients. 
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Table 3: Mean Plate Scale Errors 

Epoch Plate scale error, e 

: 1996 -0.0068 ± 0.000l3% 

1997 +0.0155 ±0.00005% 

1998 -0.0730 ±0.00018% 

1999 -0.2701 ± 0.00022% 



We have verified, using our simulations, that for 
plate scale errors of a fraction of a percent Equa- 
tions 25 and 27 are satisfied. These simulations 
show that the plate scale induced frequency error 
scales with e according to Equation 25 to better 
than 0.01% on average, with a scatter of 0.3%, for 
plate scale errors up to 0.12%. Simulations with 
e = 0.16% show a departure from Equation 25 at 
the 2% level for I > 700, while simulations for 
e = 0.27% show a departure from Equation 25 at 
the 5% level for I > 450, as illustrated in Fig. 16. 

Figure 17 illustrates how well Equation 25 is 
satisfied for the observed frequency differences. It 
shows that for the modes whose frequency is above 
1.9 mHz (i.e., n > 2) the relation in Equation 25 
is well satisfied, hence the change in frequency dif- 
ferences between epochs can be explained by the 
change in plate scale. 

On the other hand our simulations show that 
Equation 27 is only satisfied to 2% on average, 
but the plate scale induced changes in a\ are 
themselves much smaller than the corresponding- 
changes in frequency. Our simulations show that 
the relation breaks down at large degree for large 
plate scale errors, as illustrated in Fig. 18. Finally, 
Figure 19 depicts the correlation of the change in 
d\ with epoch, and shows that in practice the scat- 
ter in the data makes it hard to establish whether 
Equation 27 is actually satisfied by the measure- 
ments. 

5.2. Effect of Image Distortion 

So far we have ignored the distortion of the ob- 
served solar image. Further progress requires the 
inclusion of the image distortion and the instru- 
ment MTF, adding one more level of complexity, 
but an element that is expected to be somewhat 
stable with time. 

A full ray tracing of the optical design of the 
MDI instrument, in the Dynamics Program con- 



figuration, predicts distortions (Scherrer, private 
communication) that can be modeled as: 

e r = — =o r ((— ) 2 -l) (28) 
r r m 

where a r = 1.1 x 10~ 3 , r is the distance from the 
detector center, Ar the image distortion, and r m 
is the observed image mean radius. 

However, this expression does not account for 
the elliptical shape of the MDI image, an addi- 
tional distortion that is also on the order of 0.1%. 
A possible explanation is that the CCD detector is 
tilted with respect to the focal plane. The shape 
of the solar image resulting from such a tilt is in- 
deed elliptical, as shown in Appendix A, and a tilt 
of ~ 2° can account for the observed 0.1% effect. 

We have computed leakage matrices that in- 
clude some image distortion. These were used to 
generate simulated power ridges that were then 
fitted as described earlier. In one set of leakage 
matrix computation we only included an image 
distortion as described in Equation 28 (hereafter 
Model 1). In the second set we included the ellip- 
tical distortion resulting from a tilt of the detector 
as well the radial distortion (hereafter Model 2). 
The resulting frequency differences, as well as the 
observed ones for 1996, are shown in Figure 20. 
The distortion of the eigenvalues by differential 
rotation was included in these calculations. 

The effect of an image distortion on the fre- 
quency differences, although similar, is different 
from the effect of a plate scale error. Indeed, the 
relation expressed in Equation 25 is not satisfied 
for the simulations that include an image distor- 
tion, as shown in the lower panel of Figure 20. 

The effect of an image distortion on the odd 
splitting coefficients is comparable in amplitude 
to the effect of a similar plate scale error, as il- 
lustrated in Figure 21. The effect on the even 
coefficients is negligible. The mean change in the 
splitting coefficient offsets is listed in Table 4. 
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Tabic 4: Mean Change in Splitting Coefficient Offsets Due to Image Distortion 



Model 1 — Reference Model 2 — Reference 



5(Aai) 


-0.3457 ±0.0686 


-0.2162 ±0.0459 


<5(Aa 2 ) 


0.0017 ±0.0019 


0.0010 ±0.0012 


S(Aa 3 ) 


0.0205 ± 0.0401 


0.0127 ±0.0290 


6(Aa 4 ) 


-0.0017 ±0.0018 


-0.0010 ±0.0011 


6(Aa 5 ) 


0.0075 ±0.0119 


0.0051 ±0.0107 


5(Aa e ) 


0.0001 ± 0.0005 


0.0001 ±0.0004 



Model 1: radial term only 



Model 2: radial term and tilt of CCD 



5.3. Effect of the Instrumental Point Spread 
Function 

Another point to consider in the leakage matrix 
calculation is the effect of the instrumental point 
spread function (PSF). Not only will the smearing 
of the image by the PSF affect the leakage matrix 
(especially at high degrees) , but it will change with 
time as the focus changes. 

To estimate the PSF of the MDI instrument and 
its variation with epoch, we ran the HGE0M proce- 
dure on MDI full-disk images. HGEDM, a procedure 
part of the GONG reduction and analysis software 
package (GRASP), returns an estimate of the az- 
imuthally averaged modulation transfer function 
(MTF, i.e., the Fourier transform of the PSF), 
following the methodology described in Toner & 
Jcfferies (1993). This method reduces the 2D solar 
image to a radial profile (by computing the mean 
over concentric annuli) and computes the Hankel 
transform of this profile. By exploiting the zero- 
crossing properties of the Hankel transform one 
recovers the true image dimensions and estimates 
the MTF. 

Figure 22 (top panel) shows the PSF for five 
different full-disk images taken by MDI on five 
consecutive years. Although their shapes are very 
similar, the width of the PSF is increasing with 
the amount of defocus, as expected (Figure 22, 
lower panel). While the PSF of MDI displays a 
substantial tail beyond 2 pixels, its core can be 
approximated by a Gaussian, with a half width 
at half maximum of 0.8 ± 0.1 pixels. The image 
sizes estimated by HGE0M are consistent with the 
the plate scale based on the MDI standard reduc- 
tion pipe-line to within .04 pixels or .008%. This 
is, of course, to be expected, since the level 1.8 



MDI data processing uses the same algorithm for 
the limb definition as the GONG procedure. 

Simulated ridge frequencies resulting from in- 
cluding the PSF in the computation of the leakage 
matrix were computed. We used a Gaussian PSF, 
a Lorentzian PSF as well as the PSF estimated by 
HGE0M. The resulting changes in the frequency off- 
sets and splitting coefficient offsets are shown in 
Figure 23 as a function of frequency and degree. 
The changes in frequency offsets are on the order 
of 0.1 ^Hz (i.e., a 5 % effect), while changes in 
Acii are on the order of 0.1 nHz (or a 2% effect). 

Variations of the width of the PSF will affect 
predominantly the amplitude of the leaks with a 
small effect on the leakage asymmetry (as shown 
in Rabello-Soares et al. 2001). The change of the 
limit power spectrum computed for 1996 and 1997 
(as shown in Fig. 24 - top panel) is qualitatively 
similar to the variation of simulated power spec- 
tra synthesized using leakage matrix computations 
that includes the effect of the PSF for two different 
widths (Fig. 24 - lower panel). Our simulations 
seem to reproduce the observed changes in ampli- 
tude, but for a change in PSF width nearly 5 times 
larger than the change observed in the correspond- 
ing azimuthal averages. This might result from the 
fact that the true PSF is not azimuthally symmet- 
ric and that the azimuthal averages we used so far 
arc in fact poor approximations. 

Schou & Bogart (1998) also noticed that the 
power spectrum of MDI images varies with az- 
imuth (see Figure 1 of their paper). Figure 25 
shows power as a function of the azimuthal an- 
gle for horizontal wavenumber kf t = 230 pixel -1 
for test images taken in 1999 at four focus set- 
tings. This figure shows not only that the power 
is a function of azimuth, but, as first noticed by 
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Duvall (private communication), the phase of this 
angular dependence changes with focus position, 
changing by nearly 90° between focus positions 4 
and 5, corresponding to a change from positive to 
negative defocus (the best focus at disk center in 
1999 was at « 4.1, as shown in the upper panel of 
Figure 12). Such a change can be explained by the 
combination of an astigmatism with a positive or 
negative defocus. Such configuration would pro- 
duce a PSF that is a function of the azimuthal 
angle (Bush, private communication). 

Tarbell et al. (1997), using phase-diversity anal- 
ysis, estimated the PSF of the MDI instrument in 
focus and in its high resolution configuration (i.e., 
when the image is magnified by a factor 3.2). This 
PSF shows a clear angular dependence, with a sec- 
ondary peak with an amplitude of 0.15, attributed 
to the residual optical aberrations of the instru- 
ment components. 

Finally, we should also point out that the best 
focus changes across the CCD by nearly one focus 
step. This is caused by some residual field cur- 
vature and will produce a radial variation of the 
PSF. The optical characteristics of the MDI in- 
strument are such that the instrumental PSF has 
to be a function of both radius and azimuth. 

5.4. Effects of Image Orientation, P and 
B Q Angles 

The spherical harmonic decomposition is in 
principle computed with respect to an orientation 
aligned with the solar rotation axis. The location 
of this axis is determined by two angles, the posi- 
tion angle, P, and the hcliographic latitude of the 
central point of the disk, B . The roll angle of the 
spacecraft is maintained to keep MDI P c g to be 
zero (where P c ff is the position angle of the north- 
ern extremity of the rotation axis with respect to 
MDI CCD detector column orientation). However, 
Cliff Toner (private communication) estimated the 
value of P e ff by inter-comparing MDI and GONG 
images, drift scans, and observations taken during 
the November 1999 transit of Mercury in front of 
the Sun. He found that P cS = 0.19° ± 0.04° us- 
ing May 2000 observation and P cS = 0.20° ±0.05° 
based on the transit of Mercury. 

The precision to which B Q is known is directly 
related to the precision of the Carrington ele- 
ments, ic and flc (i.e., the inclination of the so- 



lar equator to the ecliptic and the longitude of the 
ascending node of the solar equator on the eclip- 
tic plane respectively). Giles (1999) showed that 
these are not perfect and that Aic = —0.091° ± 
0.012° and ACt c = -0.1° ± 0.1°. The error on B Q 
is therefore no larger than 0.1°. 

To assess the effect of an error in the image ori- 
entation used for the spherical harmonic decom- 
position, we computed simulated frequency differ- 
ences by fitting simulated ridges based on a leak- 
age matrix computed with either an error in P e g 
or B D of 0.25°. Figure 26 compares the observed 
frequency differences for 1996 to the theoretical 
differences resulting from including such orienta- 
tion error. The effect of a quarter of a degree error 
is unnoticeable for the frequency differences (less 
than 1% of relative change of the offsets, for either 
angle) , while an error in P e s affects noticeably the 
odd splitting coefficients. This should not come 
as a surprise: an error in the orientation of the 
image with respect to the rotation axis will affect 
the relative mode power as a functions of rajt - 
in a way conceptually similar to the distortion of 
the eigenf unctions by the differential rotation, see 
Section 3.2 - and will therefore affect the ridge 
rotational splittings. The smooth dependence of 
Aai with I for a conservative 0.25° error in P e g 
seen in Figure 26 is not readily seen in the MDI 
measurements. But this effect is small and is most 
likely masked by the random scatter and other sys- 
tematics present in the measurements. 

6. Results 

6.1. Error Budget 

Table 5 compares the range of uncertainties re- 
sulting from ridge fitting, i.e. the formal error bars 
of the ridge fitting, to the magnitude of the ridge 
to mode corrections, for both the frequency and 
the rotational splitting coefficients. While the un- 
certainties resulting from the fitting scale with the 
ridge width, Table 5 shows that the required cor- 
rections are about an order of magnitude larger 
than the uncertainties. On the other hand, to pro- 
duce useful high-degree mode parameters we need 
to be able to compute these corrections with a pre- 
cision substantially smaller than the uncertainties. 

Table 6 presents an error budget for the ridge 
to mode correction, breaking it up as a function 
of the different elements used to compute the cor- 
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Table 5: Comparison between Uncertainty and Ridge to Mode Correction 





frequency 


frequency splittings [nHz] 


a\ a 3 a 5 


uncertainty^ 
largest correction^) 


0.1 - 0.6 
4.0 


0.3 - 1.0 0.2 - 0.8 0.15 - 0.7 
10.0 5.0 2.0 



( a ) given range corresponds to a ridge width range of 10 - 60 ^iHz 
W for the 100 < £ < 900 range 



rection. The listed typical values are based on our 
best estimates of these different elements resulting 
from the current analysis. 

Now that these effects are better understood 
and that we better understand the MDI instru- 
ment imaging imperfections, we can - and will in 
the near future - reprocess the Dynamics data, 
and include the correct image scale and the known 
image distortion in the spatial decomposition. Ta- 
ble 7 presents an error budget projected for the 
best achievable values of these elements. The 
sobering conclusion of this exercise is that our es- 
timate of the achievable precision of the frequency 
and splitting coefficients corrections are still com- 
parable to the fitting uncertainties. 

For the frequency corrections, the dominant 
contributions come from the plate scale error and 
the image distortion. If both could be known to 
an accuracy of 10~ 5 , the precision of the correc- 
tion would be further reduced to 0.06 /zHz, a use- 
ful value since it would then become substantially 
smaller than the uncertainties. 

The precision of the fit of an ellipse to the limb 
of the image is likely to be 10~ 5 (i.e., fitting 5 
parameters to some 3000 limb positions). But it 
would be presumptuous to assert that the algorith- 
mic definition of the limb can be accurate to that 
level. Characterizing the instrument image distor- 
tion to the 10 -5 level is also very challenging. The 
image distortion in the MDI instrument is only 
at the 10~ 3 level (corresponding to displacements 
smaller than 0.95 pixels). We thus need to know 
the distortion with a relative precision of 1% (i.e., 
1/100 of a pixel) to reach the 10~ 5 level. 

For the splitting coefficient a\ corrections, the 
dominant contributions come from the P angle, 
and might prove to be a fundamental limitation for 
the use of high degree rotational splittings. If we 



could reduce the error in the position angle by an 
additional order of magnitude (i.e., to 0.005°), the 
precision of the correction in a\ would be further 
reduced to less than 0.1 nHz, a value now smaller 
than the fitting uncertainties. At that level, con- 
tributions from the precision of the plate scale, 
differential rotation and image distortion to this 
correction are comparable to that from the image 
orientation. 

The determination of the shape of the differ- 
ential rotation just below the surface can be im- 
proved by including high degree splittings, but 
that shape must be known to correct these split- 
tings, as they are the result of ridge fitting. It 
is likely that a bootstrapping method — or some 
other optimized methodology — can be devised to 
disentangle this problem. This still leaves the un- 
certainty of the image orientation as a substantial 
contribution. 

On the other hand, if we focus on the preci- 
sion of relative changes, i.e. differences between 
epochs, we can be more optimistic. Indeed, the 
residual systematic errors introduced by the im- 
perfect corrections are likely to be nearly constant 
with time. 

6.2. Changes with Solar Cycle 

The four Dynamics data epochs analyzed for 
this study were taken during the rising part of so- 
lar cycle 23. The first set, observed in 1996 cor- 
responds to a period when the sun's magnetic ac- 
tivity was nearly at its minimum, while for the 
forth data set, taken in 1999, the sun was near 
maximum activity. The average relative sunspot 
index for the four epochs is shown in Table 8, i.e. 
the ratio of the sunspot number to the maximum 
value of cycle 23, reached in 2000. We can thus 
study the changes of high-degree p-mode parame- 
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Table 6: Error Budget for the Ridge to Mode Correction 





typical 


corresponding change in 


parameter 


values 


frequency [yuHz] 


oi [nHz] 


plate scale 


0.01 - 0.27 % 


0.2 - 6.0 


0.06 - 1.6 


differential rotation 


4.4 - 6.2 % 


0.01 - 0.8 


0.4 - 0.7 


image distortion 


0.1% 


1.0 - 1.5 


0.2 - 0.4 


instrument PSF 


(*) 


0.05 - 0.15 


0.05 - 0.1 


image orientation, P e ff 


0.25° 


0.03 


0.01 - 2.0 


image orientation, B Q 


0.1° 


0.004 


0.012 



(*) Changes resulting from using a Gaussian versus a Lorentzian PSF 



Note also that the values tabulated are unsigned magnitude 



Table 7: Optimal Error Budget for the Ridge to Mode Correction 





achievable 


change in 


parameter 


values 


frequency [/iHz] 


a\ [nHz] 


plate scale 


0.01 % 


0.220 


0.06 


differential rotation (*) 


0.03 % 


0.054 


0.045 


image distortion 


0.01 % 


0.150 


0.04 


instrument PSF 


(t) 


0.015 


0.01 


image orientation, P e g 


0.05° 


0.006 


0.40 


image orientation, B Q 


0.02° 


0.001 


0.002 


correction estimated precision 




0.272 


0.409 


fitting uncertainties 




0.1 - 0.6 


0.3 - 1.0 



(*) values scaled from our worst case, see discussion in text 



(f) l/10th of corresponding values in table 6 



ters with the solar cycle by comparing results from 
our four data sets 

6.2.1. Frequency 

To estimate high-degree mode frequency changes 
with the solar cycle, we have corrected the ridge 
frequencies for the plate scale error (as described 
in Section 5.1, Eq. 23). We thus assume that 
all other instrumental effects are either constant 
with time or that the impact of their variation 
is small, effects such as the image distortion, the 
instrumental PSF, and the error in the P and B Q 
angles. 

Figure 27 shows frequency changes with respect 
to 1996, as a function of frequency and degree. 
These plots show that these differences increase 
with frequency, degree and sunspot number. To 
compare these changes to other investigators mea- 



surements, we show in Figure 28 mean frequency 
changes around 3 mHz (i.e., the frequency change 
averaged over the 2.6 to 3.4 mHz frequency range), 
with respect to 1996 values, as a function of de- 
gree. These are compared to estimates based on 
intermediate-degree observations (£ < 150) taken 
during the raising part of cycle 22 by Libbrecht 
& Woodard (1990), a period that corresponds to 
a change in the sunspot number index of 85, and 
the raising part of cycle 23 by Howe, Komm & Hill 
(1999), a period that corresponds to a change in 
the sunspot number index of 56. A linear extrap- 
olation with respect to degree of the shifts based 
on intermediate-degree observations agrees rather 
well up to i = 600. This confirms that these shifts 
scale with the inverse of the mode mass, indicat- 
ing that the source of the perturbation responsible 
for these shifts is located very close to the surface. 
There are also claims of a sharp drop of the fre- 
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Table 8: Relative Average Sunspot Index 

Sunspot number change 



Year 


Relative Index 


with respect to 1996 


1996 


7% 




1997 


16 % 


11 


1998 


46 % 


47 


1999 


72 % 


79 



quency shift at high frequencies. In Libbrecht & 
Woodard (1990), the drop is seen around 4 mHz, 
while in our measurements it seems to occur at 
higher frequency (i.e., around 4.5 mHz; see top 
panel of Figure 27) . One also notices that besides 
a linear increase of the frequency shift with degree, 
there is a bump centered on ^=700. 

We have checked our assumption that all other 
instrumental effects can be neglected for our esti- 
mate of the frequency shifts shown in Figure 27. 
Figure 29 presents residual frequency offsets with 
respect to 1996, namely the the quantity Az/(Pj) — 
Ai^(1996), where Azv = v' — v, v' is the ridge fre- 
quency corrected for the plate scale error, and 
v the mode frequency computed from resolved 
modes at low and intermediate degrees. This fig- 
ure shows that while there remains a residual dis- 
crepancy, there are no systematic changes with 
epoch, except for the values at low frequencies 
[y < 1.9 mHz). We can thus be confident that the 
variations seen Figure 27 are dominated by solar 
cycle effects, except for the increase in the fre- 
quency shift at low frequency, that is most likely 
an artifact resulting from some residual systematic 
errors. 

6.2.2. Splitting Coefficients 

The splitting coefficients were corrected for the 
plate scale error and their change with solar cy- 
cle estimated. Figure 30 shows the changes of the 
odd coefficients (a\ and af) with respect to 1996, 
for different epochs as a function of frequency and 
degree. The plots show that the changes of the a\ 
coefficient increase with frequency, while display- 
ing only a weak dependence on degree. Although 
the shift increases with activity, it is larger in 1998 
than in 1999. The changes for a% are very small, 
and the scatter large enough to be consistent with 
zero. 



To check if instrumental effects, other than the 
plate scale error, can be neglected and arc not 
causing the changes seen in a\, we again com- 
puted the residual ai offsets with respect to 1996, 
namely the the quantity Aai(Pi) — Aai(1996), 
where Aai = a'i — ai, a[ is the ridge fitting value 
corrected for the plate scale error, and a\ the mode 
fitting value. While there remains a residual dis- 
crepancy, there appears to be no systematic vari- 
ations with frequency, degree, or epoch. 

Similarly we have looked at the changes in the 
even splitting coefficients and in the residual off- 
sets with respect to 1996, also as a function of fre- 
quency and degree. The changes in a-i appear sys- 
tematic with frequency, degree and epoch, while 
the scatter in the residual ai offsets are large and 
inconclusive. 

In Figure 31, the mean splitting coefficient 
change and the mean residual offset, correspond- 
ing to the average over all available modes, are 
plotted as a function of epoch. Since the mean 
splitting coefficient changes track the mean resid- 
ual offsets (except for a%), we must conclude 
that the changes with epoch result from resid- 
ual systematic errors, and are not an indication of 
changes due to solar activity. 

6.2.3. Asymmetry Parameter 

Unlike the other mode parameters, the estimate 
of the asymmetry parameter is not much affected 
by the fact that we observe and fit ridges instead 
of individual modes. Korzcnnik (1999) pointed 
out that the asymmetry of the ridge seems to be 
the same as the asymmetry of the modes that are 
blended into that ridge. We found from our simu- 
lations that these differences are smaller than 10% 
up to I — 700 when including the various instru- 
mental effects described in Section 5, as illustrated 
in Figure 32. 
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Figure 33 shows the measured asymmetry coef- 
ficients for the 1996 Dynamics observations. The 
asymmetry is mainly a function of frequency, ex- 
cept for the low order modes (n = and n = 1) 
whose asymmetry is larger. In Figure 34 we com- 
pare our determinations of the mode asymmetry 2 
with values obtained for low-degrees (t < 3) using 
GOLF data (Thiery et al. 2000) and for interme- 
diate degrees using GONG observations (Basu & 
Antia 2000). Our measurements agree very well 
with the GONG determinations in the 2.4 to 3.8 
mHz range. 

In Figure 35 we compare the asymmetry for dif- 
ferent epochs, shown as a function of frequency 
and degree. Except for the /-mode, there seems 
to be no change with epoch, while the asymmetry 
of the /-mode appears to have changed systemat- 
ically, suggestive of a solar activity effect. 

6.2.4- Linewidth 

Measurements of the mode linewidths provide 
clues to the damping mechanisms of the solar os- 
cillations and are directly related to the lifetimes 
of the modes. Previous measurements have shown 
that the linewidths increase both with frequency 
and degree (Hill et al. 1996). 

Korzennik (1990) proposed a simple formula to 
estimate the mode linewidth from the the ridge 
linewidth, namely: 

Tddgo — Tmodo + ^WF + ( s ^£ ) 2 (29) 

where Twf is the width of the window function, 
sdvjdi corresponds to the effective leakage ma- 
trix width in frequency space, where the empirical 
coefficient s ~ 0.5. 

Figure 36 shows the estimated mode linewidth 
for the 1996 Dynamics observations, as well as 
linewidths estimated for low and intermediate 
degrees (0 < i < 250) based on 360-day long 
time-series of MDI Structure Program (Schou 
1999). Preliminary estimates of linewidths for 
high-^ modes (700 < I < 1800) based on high- 
resolution MDI observations (Duvall et al. 1998) 
are also included. Note that modes whose cor- 
rected linewidth are smaller than the width of the 



2 we plotted the asymmetry parameter B, as defined in 
Nigam & Kosovichev (1998), where B = , , 



window function are not shown in Figure 36, since 
the correction scheme of Eq. 29 is not valid. 

Our simple minded linewidth correction scheme 
connects reasonably well the MDI Structure 
linewidths (i.e., low and intermediate mode fit- 
ting) with the corrected MDI Dynamics linewidths, 
but we could not reconcile the Duvall et al. (1998) 
determinations based on high-resolution MDI ob- 
servations. We do see in our measurements the 
trend seen in Duvall et al. (1998) at higher de- 
grees, namely that the linewidth of the /-modes 
increases with degree faster than the linewidth 
of the pi modes. Our measurements suggest a 
cross over around I = 1000, while the Duvall et al. 
(1998) estimates imply a cross over below £ = 700. 

Figure 37 shows the relative change of the cor- 
rected linewidths with epoch - using 1996 as ref- 
erence - as a function of frequency. The magni- 
tude of these changes - on the order of 5 to 10% 
- do not correlate with solar activity, nor do they 
scale with the plate scale error. The fact that the 
linewidths of all three later epochs have changed in 
a systematically similar way with respect to 1996 
is intriguing and suggestive of residual systematic 
errors at the 5 to 10% level. 

7. Conclusion 

We believe that we have shown that one can 
construct a physically motivated model (rather 
than some ad hoc correction scheme) of the ridge 
power distribution that results in a methodology 
that can produce an unbiased determination of 
high-degree modes, once the instrumental charac- 
teristics are well understood and precisely mea- 
sured. 

We have gained substantial insight in the un- 
derstanding and modeling of the instrumental ef- 
fects of the Michelson Doppler Imager, including 
plate scale error, image distortion, point spread 
function and image orientation. Now that we bet- 
ter understand the MDI instrument imaging im- 
perfections, we can - and will in the near future - 
reprocess the Dynamics data, and include the cor- 
rect image scale and the known image distortion 
in the spatial decomposition. 

We have produced an error budget for the char- 
acterization of high degree mode parameters. Be- 
yond MDI Dynamics observations, this has di- 
rect applications on the characterization of the 
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GONG+ instruments and the reduction of the 
GONG+ observations as well as on the character- 
istics of the Helioseismic/Magnetic Imager (HMI) 
instrument planned for the Solar Dynamics Ob- 
servatory (SDO) mission, if one wishes to make 
effective use of the high degree modes accessible 
to these experiments. 

Last but not least, we also see changes with 
solar activity in the high degree frequencies and 
the asymmetry of the /-mode. Our attempt 
to see changes in splitting coefficients and mode 
linewidths remained inconclusive. 
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A. Elliptical Image Distortion 

The observed elliptical shape of the solar image is likely to be the result of an unwanted tilt of the detector 
plane with respect to the focal plane. This appendix shows how the combination of the distortion resulting 
from this tilt with the optical package distortion (or cubic distortion) can be estimated. First the cubic 
distortion is applied to the undistorted image (x ,yo) 

x = x + e r x (Al) 
y = yo + € r ya (A2) 

where e r — e r (r ) is given by Equation 28 and r = \J x^ + y%, while (x,y) are the reference coordinates in 
the focal plane that would correspond to the coordinates for a perfectly aligned detector. 

To model the image distortion due to a tilt of the CCD, the reference coordinate system (x, y) is first 
rotated around the optical axis (the z-axis) by an angle (3, which defines the direction around which the 
detector is tilted (the ir'-axis): 

x' = x cos (3 + y sin (3 (A3) 
y' = — xsin(3 + y cos/3. (A4) 

The effect of a tilt by an angle a around the x'-axis is then computed using a simple pinhole camera model. 
It is easy to show that such a model leads to the following equations: 

x" = x'{l + ^a)-x'^- (A5) 
/ 4 

V" = y'(l + ^ + V ja)-y'^- r -fa (A6) 

where / is the effective focal length and r m the observed image mean radius, both of which are measured 
in the same units as x' and y' . Note that the second term in both equations assures that the distorted and 
undistorted images have the same mean radius, while the third term in the second equation keeps the center 
of both images unchanged. 

Finally, the coordinate system is rotated back, around the z-axis, to be co-aligned with the original 
reference coordinates system: 

x'" = x"cos/3-y"sin/3 (A7) 
y'" = x" sin /3 + y" cos /3 (A8) 

and the image distortion is given by: 

S x = x'" - x (A9) 
S y = y'" - y (A10) 

where (x"',y"') are the observed (i.e. actual) coordinates and (x ,yo) are the undistorted ones. 

Due to the tele-photo nature of the MDI optical package (Scherrer et al. 1995), the effective focal length 
/ defined in Equations A5 and A6 is not the one that reproduces the image size. A complete ray trace of 
the MDI optical system (Scherrer, private communication) produces an estimate for / of 13000 pixels. 

Various observational methods have been used to determine the parameters of the distortion, some of 
which are described below. Results are summarized in Table 9. With (3 = 56° and a = 2.59°, the model 
reproduces rather well the nearly elliptical images of MDI. The resulting shape of the limb matches an ellipse 
to 10~ 6 (to be compared to the 0.1% departure from a circle), and the resulting image distortion pattern is 
shown in Figure 38. 
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An extensive series of ground tests were completed prior to launch to verify the performance of the MDI 
instrument (Zayer et al. 1993). In one of them, a grid target was imaged to evaluate, amongst other things, 
the image distortion. The distortion was determined by measuring the change in shape of the grid target as 
it was imaged onto different positions within the field of view. This distortion was fitted to the parameterized 
representation expressed by Equations A9 and A10, leaving a, / and a r to be freely adjustable parameters. 
The fit lead to a r some 11% higher and a some 5% lower than expected. This fit also lead to a value 
of / some 50% higher than predicted. The higher effective focal length is surprising since ray tracing is 
expected to produce a reliable value for this parameter. While these are pre-flight measurements, and thus 
cannot account for changes that could have occurred during deployment of the satellite, the nature of the 
manufacturing and assembly of the optical package make such post-launch changes highly improbable. On 
the other hand it is possible that some subtle aspect of the ground test has not been accounted for. 

In flight off-center continuum intensity images were acquired in February 2001 — with MDI operating in 
focus position 3 — to further study the image distortion. When the solar images are centered in different 
parts of the detector, the distortion will affect them in a different way and the shape of the limb will change 
with each image position. Using again the image distortion parameterized by Equations A9 and A10, the 
predicted shape of the limb can be computed as a function of image position, and compared to the observed 
limb profiles. The best match between predicted and observed limb shapes was achieved for a = 2.6°, 
13 = 56° and / = 12940. 

Finally, images obtained during rolling maneuvers of the spacecraft in March 1997 and November 2001 
were also used to characterize the distortion. These images were averaged for each roll angle and rotated 
back to P c g — 0. The resulting images, which are dominated by the supergranulation, are then used to 
estimate the differential distortion by correlating each image with an image in the center of the time range. 
In this case / was kept fixed while the same distortion model was again fitted to the data. Two estimates 
were made using these images. The first (hereafter method 1) is made using the offset of the center of the 
images, while the other (hereafter method 2) is made using the stretching of the images. The results are 
inconsistent as summarized in Table 9. In particular method 2 always produces higher values of a than 
method 2. Again, a larger value of / would be required to reconcile results from both methods. 

Table 9: Summary of the Characterization of the Model of MDI Image Distortion 



Method 


a 


P 


/ 


a r 


Match limb shape 


2.59° 


56° 


Fixed 


Fixed 


Pre-flight target grid 


2.45° 




19604.4 


1.2 x 10~ 3 


Image offset 


2.6° 


56° 


12940.0 


Fixed 


1997 roll, method 1 


1.80° 


52° 


Fixed 


Fixed 


1997 roll, method 2 


1.97° 


59° 


Fixed 


Fixed 


2001 roll, method 1 


1.71° 


53° 


Fixed 


Fixed 


2001 roll, method 2 


2.43° 


50° 


Fixed 


Fixed 



Although the parameters resulting from the different methods used to characterize the image distortion of 
MDI do not agree as well expected, the distortion pattern resulting from a tilt of the detector is a rather good 
model of the observed distortion. The fact that the model parameters are so different when using different 
methods might be attributed to the complicated instrumental PSF of the MDI optical package that is not 
fully known and therefore cannot be taken into account. In particular the different observational methods 
described above use features of the image at different spatial scales (or spatial wavenumbers) to estimate 
the image distortion. The PSF may also explain why the ellipticity of the images changes significantly when 
the images are strongly defocused. Another possible explanation may, of course, be that other types of 
distortions, such as non square pixels, are present. 
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Fig. 1— MDI sectoral power spectra, for Dynamics 1999 and for (^,n) = (100, 6), (100, 10) and (300,2), 
computed using a 9th order sine multi-taper. At low degree and low order (top panel) individual modes 
are well resolved: the isolated peak at 2970 yuHz is the target mode while the others are the spatial leaks. 
Otherwise the peaks are blended: the middle panel shows the spectrum for the same degree but at higher 
order where mode blending results from larger mode widths due to shorter life times. The bottom panel 
shows the spectrum at roughly the same frequency as in the top panel but at a higher degree, where mode 
blending results from a smaller mode spacing (dv/dt) as well as shorter life time. 
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Fig. 2. — Synthetic (dashed line) and observed (continuous line) sectoral limit spectra, Dynamics 1996, for 
n = 2 and £ = 11 and £ = 122 in the top and bottom rows respectively. Synthetic spectra computed using 
leakage matrices without and with the inclusion of the horizontal component are plotted in the left and right 
columns respectively. 
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Fig. 3. — Perturbation of the leakage matrix by the solar differential rotation: Ge/>(m), for t = 400,500 
and 600, at n = 3, see definition in the text. The panels on the left show the Gt^i coefficients as a function 
of p = (£ — £')/2 and m (dash lines and dotted lines correspond to positive and negative contour levels 
respectively), while the panels on the right show profiles at fixed m, for the m values indicated by the 
horizontal lines on the corresponding left panel. 
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Fig. 4. — Synthetic power spectrum computed for the n = 2, i = 500, m — —300 mode without and with the 
inclusion of the effect of the distortion of the eigenfunctions by the differential rotation (top and bottom panel 
respectively), and neglecting the horizontal component of the leakage matrix. The dot-dashed lines show 
the power spectra that would be observed if the mode linewidth was one tenth of the observed value. The 
vertical dash and dot lines indicate the mode frequency and the peak of the ridge respectively — displaced 
in this case by some 9.7 /zHz. 
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Fig. 6. — Differences between ridge and mode frequencies, and ridge and mode splitting coefficients, as 
a function of frequency, for the Dynamics 1996 data (black points). These differences arc compared to 
theoretical differences (for the same set of modes) computed from simulated power spectra without and with 
the inclusion of the effect of the distortion of the eigenfunctions by the differential rotation (green and red 
points respectively). 
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Fig. 7. — Plot similar to Fig. 6: differences between ridge and mode frequencies, but with theoretical 
differences extending to spherical harmonic degrees up to £ = 900. 
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Fig. 8. — Plot similar to Fig. 7: differences between ridge and mode frequencies, but shown as a function of 
spherical harmonic degree, I. 
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Fig. 9. — Comparison of the latitudinal dependence of three different rotation profiles at the solar surface, 
used to estimate the effect of the differential rotation on the distortion of the eigenfunctions. Model A 
corresponds to somewhat arbitrary values (B 2 = —75 nHz, B4 = —50 nHz); Model B to values derived from 
Schou ct al. (1998) (B 2 = —51 nHz, B4 = —81 nHz); model C to values derived from Snodgrass & Ulrich 
(1990) (B 2 = -77 nHz, B 4 = -57 nHz). 
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Fig. 10. — Changes in the ridge to mode frequency and splittings coefficients offsets caused by the distortion 
of the eigenfunctions by differential rotation, corresponding to the different rotation profiles illustrated in 
Figure 9. 
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Fig. 11. — Top panel: observed q = j3/L estimated using the Dynamics 1997 data set using S x (diamonds) 
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Fig. 12. — Top panel: variation of the MDI mean focus as a function of time, as computed by Bush, Chu, & 
Kuhn (2001). The horizontal lines indicate the focus position of the instrument at the time, and therefore 
the difference between the measured focus and the instrument focus position gives the amount of defocus 
for each Dynamics data set (vertical dashed lines). Bottom panel: resulting plate scale error, i.e., relative 
difference between the observed mean solar image radius (i? m can), calculated from each image using HGEDM, 
and the radius determined using the initial plate scale (i? SU n, see text for a more detailed description). 
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Fig. 13. — Comparison of observed and simulated sectoral limit power spectra (top and bottom panels 
respectively). The top panel compares Dynamics 1998 to 1999 spectra while the bottom panel compares 
synthetic spectra computed using leakage matrix without and with the inclusion of a plate scale error of 
—0.1% (i.e., solar image radius is 0.1% larger than the one assumed in the spatial decomposition). 
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Fig. 14. — Comparison of the differences between ridge and mode frequencies, and differences between ridge 
and mode odd splitting coefficients, as a function of frequency, for all the Dynamics data sets. 
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Fig. 15. — Theoretical differences between ridge and mode frequencies (left panel) and between ridge and 
mode a\ splitting coefficients (right panel), shown as a function of frequency, resulting from fitting synthetic 
spectra computed with leakage matrices with various plate scale errors (e = 0,0.036,0.071,0.120,0.160 and 
0.270%). 
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Fig. 16. — Relation ((ei, e 2 , €3) (see Eq. 25) that accounts for changes in the theoretical differences between 
ridge and mode frequencies resulting from different plate scale errors. The dotted lines indicate the value 
corresponding to the combination of plate scale errors given in Eq. 25. This relation is illustrated with three 
combination of plate scale errors to show that it breaks at large degrees for large plate scale error. 
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Fig. 17. — Relation £(Pi,P2,P3) (Eq. 25) for the four possible combination of Dynamics observations, as 
a function of frequency (left panel) and degree (right panel). The horizontal dash lines correspond to the 
expected values (see right-hand-side of Eq. 25). This confirms that the changes seen in Fig. 14 are primarily 
due to the change in the plate scale error. 
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Fig. 18. — Relation between Aa\ and £1,62 (see Eq. 26) that expresses the changes in the theoretical 
differences between ridge and mode a\ splitting coefficients resulting from different plate scale errors. The 
dotted lines correspond to the value predicted by Eq. 27. This relation is illustrated with three combination 
of plate scale errors to show that it breaks at large degrees for large plate scale error. 
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Fig. 19. — Average change in the observed differences between ridge and mode a\ splitting coefficients (dia- 
monds) as a function of change in plate scale errors, for all possible Dynamics observation combinations. The 
error bars correspond to the standard deviation about the mean. The crosses and stars symbols correspond 
to our simulations (all the modes or only the modes with degree up to 250 respectively) and the dotted line 
illustrates the relation expressed by Eqs 26 and 27. 
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Fig. 20. — Top panel: Comparison of observed and theoretical differences between ridge and mode frequen- 
cies. The effect of including an image distortion (in green and blue) in the leakage matrix is compared to 
our reference model (in red) and to the observations (in black). The bottom panel shows that the effect of 
the image distortion is different from an image scale error: the relation expressed by Eq 25 is not satisfied. 
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Fig. 21. — Changes in odd splitting coefficients offsets for two image distortion models with respect to our 
reference model. Model 1 includes only a radial image distortion (green points), Model 2 includes both the 
radial term and the tilt of the CCD (blue points). Horizontal lines are the respective means. 
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Fig. 22. — Top panel: estimates of MDI azimuthally averaged point spread function (PSF), computed from 
four images corresponding to each consecutive year of the Dynamics program. These are compared to the 
Airy function corresponding to the instrument aperture diffraction limit and to a Gaussian and a Lorentzian 
with a matching half width at half maximum (HWHM) . Bottom panel: HWHM of the azimuthally averaged 
point spread functions as a function of image defocus. 
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Fig. 23. — Change in the theoretical differences between ridge and mode frequency (top panels) and and 
between ridge and mode a\ splitting coefficient (bottom panels) , with respect to a reference function 
of frequency (left panels) or degree (right panels). The reference case uses a leakage matrix computation 
that does not include the effect of the instrumental PSF. This is compared to the cases where the leakage 
matrix computation includes a) a Gaussian PSF (HWHM = 0.71 pixel, black points), b) a Lorentzian PSF 
(HWHM = 1.0 pixel, red points), or c) the PSF estimated by HGE0M using a Dynamics 1996 image (green 
points). These changes are on the order of 10% for the frequency differences and 2% for the oi coefficient 
differences. 
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Fig. 24. — Observed (top panel) and synthetic (middle panel) sectoral limit power spectra for n — 2 and 
t = 92. The observed spectra correspond to the two consecutive years with the largest change in PSF width 
(see Fig 22), that might account for the observed change in leaks amplitude. The three synthetic spectra 
correspond to leakage matrices computed using no PSF, a Gaussian PSF with a HWHM of 0.71 pixel (a 
good approximation of the measured azimuthally averaged PSF) and a Gaussian PSF with a HWHM of 4.0 
pixel (a substantially wider PSF). Only when very wide PSF is included does our simulations produce a 
noticeable change in the amplitude of the leaks. 
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Fig. 25. — Observed power spectra as a function of azimuthal angle for a horizontal wave number, kh, of 
230 pixel^ 1 , computed from 1999 MDI images taken at four different focus positions, i.e. focus steps 3, 4, 5 
and 6. The best focus position at the time was determined to be at focus step 4.1 (see Fig 12). The angular 
dependence of these power spectra indicates that the instrumental PSF is not uniform with respect to the 
azimuthal orientation on the detector. Note how the non-uniformity increases with increasing defocus and 
the phase shift of the azimuthal dependence between either side of the best focus position. 
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Fig. 26. — Comparison of observed and theoretical differences between ridge and mode frequencies and 
between ridge and mode odd splitting coefficients. The effect of including an error in the solar rotation axis 
orientation in the computation of the leakage matrix (5P c s — 0.25°, in green, and 5B Q = 0.25°, in blue) is 
compared to our reference model (in red) and to the observations {Dynamics 1996, in black). 
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Fig. 27. — Binned frequency changes with respect to values for 1996, after correction for respective plate 
scale errors. The differences were averaged over bins 100 /iHz wide in frequency (top panel) and 20 degree 
wide (bottom panel); the error bars represent the standard deviation of the mean. 
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Fig. 28. — Mean frequency changes around 3 mHz (i.e., the frequency change averaged over the 2.6 to 3.4 
mHz frequency range) - with respect to 1996 values - as a function of degree. These values are compared to 
linear extrapolations of comparable determinations based on Big Bear Solar Observatory observations taken 
between 1986 and 1988 (Libbrecht & Woodard 1990) and on GONG observations taken between March 1996 
and July 1998 (Howe, Komm & Hill 1999). 
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Fig. 29. — Residual frequency offsets with respect to values for 1996, as a function of frequency and degree. 
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Fig. 30. — Changes in a\ and 03 splitting coefficients - relative to values for 1996 - for different epochs, 
as a function of frequency and degree, after correction for respective plate scale errors. The changes were 
averaged over bins 100 /iHz wide in frequency (top panel) and 20 degree wide (bottom panel); the error bars 
represent the standard deviation of the mean. 
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Fig. 31. — Mean changes in odd and even splitting coefficients with epoch (circles) compared to mean changes 
in residual offsets (stars, see text). Since both sets track each other (except for ae), we are led to conclude 
that the changes are due to systematic errors. 
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Fig. 32. — Comparison of the ridge asymmetry resulting from fitting synthesized ridges with the mode 
asymmetry used as input. The different colors represent results for different simulations. Up to I = 700 the 
ridge asymmetry is within 10% of the mode value. The large excursions are for low order modes (i.e., for 
n — and 1). Also, the feature at v = 4.7 mHz occurs where the asymmetry goes through zero. 
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Fig. 33. — Ridge asymmetry for the 1996 Dynamics epoch, as a function of frequency and degree. Error 
bars are shown only for every 5th data point, the low order modes are color coded. 
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Fig. 34. — Comparison of the profile asymmetry estimates from f 996 MDI Dynamics data (this work) with 
published values from GOLF and GONG observations. The MDI data were binned in frequency for n = 0, 1 
and n > 2 separately. The error bars for the MDI and GONG values represent the scatter of the data. 
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Fig. 35. — Changes of the asymmetry parameter - with respect to 1996 - for different epochs, as a function 
of frequency and degree. Changes for the f-modes are indicated by diamonds. 
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Fig. 36. — Corrected MDI linewidths (circles, this work) compared to MDI low and intermediate degree values 
derived from 360-day-long Structure data (crosses) and to very high degree corrected linewidth estimates 
from MDI in high-resolution mode (stars) as published by Duvall et al. (1998). Uncorrected linewidths 
are shown as dots. The horizontal lines correspond to the FWHM of the window functions for the MDI 
observations. 
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Fig. 37. — Raw and binned linewidth relative changes (dots and circles respectively), with respect to 1996, 
as a function of frequency, and restricted to modes common to all epochs. The ridge linewidths were first 
corrected according to Equation 29. 
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Fig. 38. — Map of the image distortion in the plane of the CCD. The distortion has been multiplied by a 
factor 50 to make it visible. The largest displacement inside the solar limb is 0.95 pixel. The top panel 
illustrates how the distorted departs from a circle (i.e., e c — (x/r) 2 + (y/r) 2 — 1) and how well it fits an 
ellipse (i.e., e e — (x/a) 2 + (y/b) 2 - 1). 
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